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Abstract: Metabolic configuration and adaptation under a range of abiotic stresses, 
including drought, heat, salinity, cold, and nutrient deprivation, are subjected to an intricate 
span of molecular pathways that work in parallel in order to enhance plant fitness and increase 
stress tolerance. In recent years, unprecedented advances have been made in identifying and 
linking different abiotic stresses, and the current challenge in plant molecular biology is deciphering 
how the signaling responses are integrated and transduced throughout metabolism Metabolomics 
have often played a flindamental role in elucidating the distinct and overlapping biochemical 
changes that occxxc in plants. However, a far greater understanding and appreciation of the 
complexity in plant metabolism under specific stress conditions have become apparent 
when combining metabolomics with other — omic platforms. This review focuses on recent 
advances made in understanding the global changes occiuring in plant metabolism under 
abiotic stress conditions v&'mg metabolite profiling as an integrated discovery platform. 

Keywords: abiotic stress; metabolic reconfiguration; metabolomics; plant; proteomics; 
transcriptional regulation 
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1. Introduction 

Pknt growth and development are significantly affected by varbus stress conditions. A direct 
consequence of a stress perturbatbn is an alteration in the metabolic behavior of the cell, leading to a 
cascade of molecular and bbchemical events that feciHtate a new steady state to be reached. Depending on 
the development stage of the plant or type of stress applied, the phenotypical output can be associated 
with reduced or delayed germination efficiency [1], decreased vegetative growth [2], delayed or 
arrested cell cycle activity [3,4], earlier or delayed transition to flowering [5,6], decreased formation 
and viability of reproductive organs (for review, [7]), or accelerated senescence [2,8]. Small molecules 
(metabolites) play an important role during these transitions and adaptations. Perhaps the most well 
documented changes are the accumulation of certain metabolites, such as proline (induced upon 
drought, cold, osmotic, and salinity stresses, [9,10]), soluble major and minor sugars (induced upon drought, 
light, cold, and nutrient stresses; see also Section 3.3), glycine betaine (induced upon salt, drought, and 
cold stresses [11],) and amine containing compounds such as putrescine (induced upon cold and 
oxidative stresses, [12-14]), ^ermine (induced upon cold, oxidative, and salinity stress, [13,14]), and y- 
aminobutyrb acid (GABA) (induced vpon salinity [15], anoxia [16], and carbon deprivation [17]). The 
osmotic corrpatibility within the cell is maintained by the accumulation of some of these metaboUtes 
(e.g. proline, glycine betaine, GABA, and su^rs), decreasing the entropy levels within the cell, and 
albwing for folded native tertiary proteins structures to be maintained. The accumulation of other 
metaboHtes (e.g. ascorbate, glutathione, vitamin Bi, and Be, [18-21]) also significantly reduces the 
harmful effect of reactive oxygen species (ROS) generated by abbtb stresses while ROS itself might 
act as an inportant messenger during stress responses (for review [22]). MetaboHtes can also serve as 
important allosteric regulators during abiotic stresses. AUosterb regulatbn involves the direct binding 
of a metaboHte to a protein, which modifies the interaction, localization, stability, or substrate affinity 
of the gene product. Furthermore, the contribution of metabolites (particularly phytohormones and 
sugars (Section 3.3)) as important messengers during stress signal transduction under adverse growth 
conditions has also been well documented (for example [23]). As the metabolome (metabolite 
complement of a cell) represents a snapshot of the prevailing biochemical state of a particular organ or 
tissue under investigation, and the current estimation of the chemical complexity across the plant 
lineage approximates to -200,000 chemical diverse fingerprints [24], the scope for identifying and 
elucidating more signaling networks upon abiotic stresses remain to be discovered. Moreover, 
metabolic activities respond to stress more quickly than transcriptional responses [25], suggesting that 
the elucidation of metabohte changes will increase our knowledge of how complex metabolic networks 
interact and how they are modified upon specific stresses [25]. Using metabobmics as an important 
diagnostic tool within the right discovery context, thus, provides a powerful means to gain novel 
bb logical insight into the metaboUc regulatory network. For this purpose, this review does not attempt 
to detail individual metabolite changes in response to abbtic stress conditbns; for a recent coirprehensive 
review covering these aspects we refer to [25,26]. In additbn, for analytical techniques associated with 
plant metabolite profiling, and data analyses and processing ppelines, excellent reviews have recently 
covered these [24—28]. Here, rather, we focus on the interactbn between, and wealth of infbrmatbn that 
can be uncovered in studying, metabobmbs combined with other molecular platforms during abiotic stress 
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responses. The overall result is a more thorough and intricate understanding of the global adaptation 
patterns and hierarchical reguktbn structures that prevail during abbtic stress conditions. 

2. Transcriptional Regulation of Metabolic Networks 

Key to our understanding of the central dogma of molecular biology is the linear progression of 
molecular events associated with it; initiated with the encoding and replication of the genomic 
information contained in the DN A of an organism. This gives rise to the prodiaction of messenger RNA 
(mRNA) that will allow for differential gene expression patterns to emerge during the process of 
transcription. The mRNAs are subsequently processed by an array of splicing and editing mechanisms 
before exiting the nucleus to the cytoplasm Within the cytoplasm, the mRNAs will interact with large 
ribosomal coirplexes which will translate the information contained on the mRNA molecules, leading 
to protein synthesis. Each protein sjmthesized can then catalyze a set of biochemical reactions leading 
to the dissipation of a given substrate (metabolite x) to its bbchemical (end) product (metabolite y). 
Until recently, this widely perceived interpretation of a metabolic pathway was commonly applied (and 
perhaps naively the only considerable viewpoint) when considering and interpreting the biobgical meaning 
of transcriptomic, proteomic, or metabolomic datasets when viewed in isolation. The pitfeUs associated 
with this assunption have especially become apparent when combining two or more -omic technologies in 
order to obtain a more holistic view of metabolism. For a number of exanples metabolic changes do 
follow the conventbnal transcriptional response. For instance, characterization of the core aerobic and 
anaerobic responses in rice have identified that decreased transcript abundances in sucrose cataboMsm 
or tricarboxylic acid (TCA) cycle reactions corresponds to decreased hexose equivalents and organic 
acids, respectively [29]. Consequently, a metabolic switch to fermentative metabolism is apparent. 

However, what happens if such a clear progression (correlation between transcript (or protein) and 
metabolite levels) is not eminent in the dataset? How do metabolomics (and associated bioinformatic 
tools) aid in resolving such (apparent) discordances, and how will this influence how we evalioate 
metabolite data in fixture? Perhaps, the first important aspect is that the unequivocal identification and 
robustness in the dataset(s) have to establish. Currently, metabolite profiling efforts primarily focus on 
hyphenated chromatography based technobgy coupled to mass spectrometry (gas chromatography 
mass spectrometry (GC MS) or liquid chromatography mass spectrometry (LC MS)) high throughput 
techniques to profile highly to moderately abundant metabolites. In order to achieve this, unequivocal 
metabolite identificatbn relies on the establishment and verificatbn of custom metabolite libraries that 
rely on retentbn time correction indexes (retentbn index, RI) and matching mass-to- charge (m/z) 
identifiers. While some search algorithms are using more quantitative means of matching spectra to 
specific metabolites, the overall generalizatbn still prevail to match 4-6 of the most unique m/z 
identifiers to a specific metabolite within a specific RI window. However, metabolites that share similar 
m/z identifiers, and elute in close proximity will either be misidentified using this methodology or 
would not be distinguished from metabolites sharing similar chemical properties. Improvements in 
recent MS technologies have led to the development of mass accuracy inst-umentation, which significantly 
enhance the elucidation of metabolite identities that share similar spectra information. Moreover, 
technologies such as nuclear magnetic resonance (NMR) provide the most unambiguous information 
regarding the structural properties of a smaU molecule. While the use of this instrumentation has several 
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limited capacities in terms of the sensitivity of detection, new applications such as coupling NMR and 
hyphenated technologies will significantly advance the profiling and identification of distinct metabolites, 
and albw for more robust metabolite data reporting. Irrproved peak resolution by using two-dimensional 
(or multi-dimensbnal) GCxGC time of flight mass spectrometry (TOF MS) has also been achieved to 
significantly iirprove separating power, and allow for a higher resobtbn capacity [30]. Several artifacts 
may also arise from derivatizatbn steps during sanple preparation of nonvolatile metabolites for GC 
MS analyses [31]. As an exairple, the conversion of arginine to the trimethylsilyl (TMS) derivative of 
ornithine has been previously documented [32], while the systematic detection and documentation of 
artifects arising partbularly from GC MS based experiments have been collated into in house databases 
(e.g. [33]). One approach to address such artifect formatbn is through the classbal addition of authentb 
standards or spiking experiments. The annotatbn and quantificatbn of metabolites can also be 
improved by different stable isotope labeling strategies that, in combinatbn with standardized database 
searches and elemental mass calculations, allows for reduced chemical background noise and increased 
confidence in the elemental formula annotation [34]. A highly sensitive, ultratast (15s) isotope-derivatization 
reagent, do-/d6-2, 4-dimethoxy-6-piperazin-l-yl pyrimidine (DMPP), has also been developed to specifically 
derivatize carboxylic analytes, which will reduce chemical noise [35], leading to a greater confidence 
in the unequivocal identification of these acidic compounds. Computational annotation methods to 
assist in refining metabolite annotations have also been developed (e.g. [36]), and will continue to 
improve the accuracy that metabolites are identified with 

Apart from metaboHte identification, any other apparent discrepancy in metaboHc changes also has 
to be rigorously evaluated in the biosynthetb context. Customized in-house and pubHcly available 
software packages have been developed for standabne metabolomic data and pathway enrichment 
analyses that greatly aid in statistically assessing whether metaboHte changes corresponds to whole 
pathway enrbhment analyses (e.g. [37-39]) or graph-clustering algorithms in order to understand the 
org9nizatbn of metaboHc fimctbnal modules [40]. While such analyses have not been strictly applied to 
plant specific abbtic stress metaboHte profiles to date, the integration of these applications greatly 
accelerates and creates confidence in the information generated from metabolite profiles to determine 
the metaboHc pathway activity associated with the particular perturbatba Other bioinformatb took to 
discriminate and interrogate metabolite data include both supervised and unsupervised visualization 
techniques, including batch leaming self- organizing map analysis (BL SOM) [41], orthogonal projectbn 
to latent structures discriminant analysis (OPLS DA), principal coirponent analysis (PCA) and Bayesian 
independent component analysis (ICA) in order to reduce the dimensionaHty or non-overlapping information 
in the data [42]. In addition, in a complementary approach, metabolomics and other — omic technologies 
could be performed in parallel and, for a combined metabolomics and transcriptomic approach, the output 
can be visualized; overrepresentation and other interrogative statistical analyses performed [43-46]. A 
number of studies have successfully integrated data from metabolomics to transcriptomics, and even 
enzyme activities [47-50], attesting to the value in using parallel approaches to understand plant 
metaboHsm (see also bebw). This review highlights several of these combined approaches in order to 
understand specific abiotic sfresses (or abiotic stress responses in general), and the advancement of 
knowledge made in understanding the molecular context of metaboHte changes through both targeted and 
untargeted profiling techniques. These metabolb changes are specificaUy discussed for their alternative 
roles in metaboHsm, ranging from sfructural modifiers of nucleb acids to inportant mediators or mdicators 
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of signal transduction mechanisms, concepts largely ignored when interpreting the context of metabolic 
changes during stress perturbations. 

3. Metabolic Regulation of Transcriptional Networks 

Any disagreement between transcriptomic and other -omic (proteomic and metabolomic) datasets 
are usually interpreted as an additional layer of regulation either occurring on a post-transcriptional 
or — trans lational level (although the direct comparison between technological platforms and sensitivity 
and analytical capacity of instrumentation should also be taken into account) [51]. The intriguing 
observation that certain metabolites can exert a transcriptional response preceding the inductbn or 
independent of its own biosjoithetic gene expressbn has, however, emerged in the last years [49,52]. 
In addition, the exogenous applicatbn of certain metabolites that modulate or mimic abiotic stress 
responses [52,53] also attribute to the coirplexity of interpreting metabolic changes in its strict pathway 
context. The direct influence of prevailing metabolite levels on gene and protein expression or activity 
is usually referred to as metabolic regulation, ranging fromsinple allosteric regulation to the emergence of 
direct interaction of metabolites with nuclear gene expression (NGE) during retrograde signaling 
(Sectbn 3.2). Regulation strategies are also most likely unidirectbnal as transcriptional regulatbn of 
metabolic changes will lead to the metabolic inductbn of transcriptbnal signal transduction pathways. 
In addition, overlapping abiotic stresses that act either synergistb or antagonistic in behavbr might 
exhibit phases of both. As an example, the cold acclimation response and subsequent interaction with 
carbon availability exhibit both transcriptional and metabolic regulation strategies [54]. Thus, metabolites 
can exert its eflfect(s) on multiple hierarchical levels of control, displaying much more versatile roles than 
previously anticipated. 

3.1. Metabolic DNA and RNA Structural Modifiers 

Metabolites are not iisiially considered to have a direct interaction with genomic information (apart 
from forming the backbone and building blocks of this macromolecule); however, molecular evidence 
suggests a direct metabolite mediated modulation of chromatin structures. In plants, the monophosphorylated 
isomer of the phospholipid phosphatidylinositol (Ptins), phosphatidylinositol 5-phosphate (PtIns5P), 
interacts and binds directly to the PHD domain Arabidopsis trithorax like factor, ATXl [55]. ATXl 
is an epigenetic factor that modulates chromatin structures by methylation of the N-terminal lysine 4 
(K4) of histone H3 (H3Kme3) [56] but acquire target specificity as only specific genes such as the 
transcription fector, WRKY70 [56], or the flowering repressor gene, FLOWERING LOCUS C [57], are 
transcriptionally activated by ATXl. In response to abiotic stresses (in particular salt, osmotic, and 
drought stress [55,58]), PtIns5P levels increase significantly, leading to the direct binding of C5^osolic 
PtInsSP to ATXl. This interactbn prevents the nuclear import of ATXl, leading to decreased 
WRKY70 expression [59]. WRKY70 is implicated as a central convergence node between jasmonic 
and salicylic acid mediated signaling networks [60], and ne^tive regulator of leaf senescence [61]. 
Such metabolite mediated control over epigenetic chromatin remodeling potentially has fer reaching 
consequences, as it can be directly linked to phenotypical outputs associated with senescence and 
flowering aberratbns during stress conditions. Moreover, small molecule moduktbn of chromatin 
structure and its subsequent effects can be experimentally empbyed to distinguish and clarify between 
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the contentious sijbject of a transient DNA methylation state (influencing seed filling and subsequent 
germinatbn efficiency) and trans generatbnal heritable epialleles (linked to an increase in recombinatbn 
fi-equency [62]), associated with the "memory" or "priming" effect that is fi-equently observed in the 
offspring of parents that are exposed to abbtic stresses. 

Another unexpected interaction for small molecules in plant metabolism involves the existence of 
metaboUte binding riboswitches. Riboswitches are mRNA elements that control gene regulation or 
resulting protein translation via a conformational change in the RNA structure ipon direct binding of 
small interacting molecules. In bacteria and archaea, more than a hundred putative metaboHte/small 
RNA (siRNA) responsive RNA sensors are predbted to exist [63], with 13 experimental^ vaHdated 
elements that bind metaboHtes specifically (namely purines (adenine, guanine), S-adenosyl-homocysteine 
(SAH), S-adenosyl methionine (SAM), tetrahydro folate, pre-queuosonei, lysine, glycine, glutamine, 
glucosamine- 6-phosphate, flavin mononucleotide (FMN), cyclic di-GMP, adenosyl-cobalamin/vitamin 
Bi2, and thiamin pyrophosphate (TPP)). In plants and fLingi, however, only the TPP riboswitch has 
been identified to date, with the other known riboswitches not conserved in the plant lineage [64,65]. The 
plant riboswitch is located at the 3' untranslated region (UTR) of the thiamin monophosphate (TMP) 
bbsynthetic gene, THIAMIN C SYNTHASE (THIQ. Jn Ambidopsis , two other enzymes, namely THl and 
THIl are also involved in TMP biosynthesis; however, in contrast to cytosolic localized THIC, occur in the 
chloroplast. TMP is subsequently dephosphorylated into thiamin, which is then pyrophosphorylated into 
TPP by cytosolic thiamin pyrophosphokinases (TPKs), TPKl, and TPK2 ([66], and references 
therein). In the nucleus, direct TPP binding to the riboswitch leads to the intron spHcing of the catalytic 
domain of THIC for possible nonsense mediated decay, leaving only the riboswitch structure and 
thereby controlling endogenous thiamin and TPP levels [66]. The mechanism of plant TPP riboswitch 
control over metaboUsm has not been shown to have a direct role in abiotic stress sensitivity. However, 
exogenous application of thiamin and TPP confer enhance stress tolerance to high light, cold, 
osmotic, salinity and oxidative stress conditions [19,67]. It is also able to rescue the ROS sensitive 
ascorbate peroxidasel mutant from oxidative stress by decreasing protein carbonylation events and 
H2O2 production [19]. Abbtb stresses increase the endogenous levels of TMP, thiamin and TPP as a direct 
consequence ofthe induction of its own biosynthetb gene {THIC, THII, THII, and TPATi) expression [19]. 
Furthermore, TPP is an essential cofactor for the tricarboxylic acid (TCA) cycle enzyme complexes, 
pyruvate dehydrogenase (PDH) and 2-oxoglutarate dehydrogenase (2-OGDH), and for the a-ketose 
transketolase (TK) homodimer of the pentose phosphate pathway. Consequently, TPP availability 
controls the rate of carbohydrate oxidation through these pathways [66], both pathways (as well as 
carbon statuses, see Section 3.3) significantly affected by abiotic stresses. While THIC promoter 
activity is controlled by the circadian clock [66], the expression ofTHI and THII are dependent on the 
stress inducible phytohormone, abscisic acid (ABA) [67]. In contrast, TPK expression shows no 
variation in response to abiotic stresses, ABA or circadian regulation [68]. In light of these observations, 
the dynamic relationship between the subcellular thiamin and TPP biosynthetic pathways and 
subsequent involvement ofthe TPP riboswitch (and the mutated riboswitch variant [66]) ipon abiotic 
stress conditbns would provide an interesting opportunity to explore the mechanistic basis for thiamin 
metabolism ipon (partbularly) oxidative stresses. The prospect of also engineering plants with customized 
riboswitches to transiently respond to changes in metaboUte levels is a novel avenue to e?q)lore in 
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future for enhancing abiotic stress tolerance without growth penalties or aberrant phenotypic variations 
usually observed for constitutive gene alterations grown under non stressed growth conditions. 

3.2. Retrograde Signaling 

Communication between organelles and nucleus is imperative to ensure a coordinated orchestration of 
avoidance or tolerance mechanisms in response to various stress stimuli. Plant organelles produce multiple 
retrograde signals to alter and coordinate nuclear gene expression (NGE) with organellar gene expression 
in order to optimize and activate defense, avoidance, or tolerance mechanisms during adverse growth 
conditbns [69-71]. Inplastids, certain metabolites, including tetrapyrrole Mg-protoporphyrin (Mg-Proto IX, 
induced during high light stress, [72]), 3'-phosphoadenosine 5 '-phosphate (PAP, induced during drought or 
hi^ light stress, [73]), methyferythritol cycbdiphosphate (MecPP, induced during wounding or hi^ Ught 
stress, [74]), ROS (induced during oxidative or high light stress, [75]), ABA, or P-cycbcitral (induced 
during oxidative or high light stress, [76]) have all been reported to accumulate in response to varbus 
abbtic stresses and inp heated as inportant signals during organellar retrograde signaling. 

The best documented plastid retrograde system is that involving intermediates of pigment sjoithesis 
identified using a forward genetic screen, identifying genome uncoupled (gun) mutants [72] in which 
the expressbn of light harvesting chlorophyll a/b binding (LHCB) is uncoipled from the functbnal 
state of chlorop lasts. Five of the six gun mutants, gun2-6, have impaired flux through the tetrapyrrole 
bbs5aithesis pathway, accumulating lower levels of Mg-Proto IX, an intermediate in chlorophyll 
bbsjoithesis [72]. Whilst initial reports established a direct correlation between Mg-Proto IX levels 
and transcript abundances, subsequent stvidies with more sensitive technologies have revealed that Mg- 
Proto IX transiently accumulate [77] while mutants or chemical inhibition of the chlorophyll 
biosynthesis pathway lack a clear correlation between nuclear gene expression and the tetrapyrrole 
intermediate [78,79]. Studies with the gain of function gwn 6-lD mutant has fLirther implicated heme 
levels, rather than Mg-Proto IX, as the most promising candidate to coordinate photosynthesis 
associated nuclear gene (PhANG) expression with chlorop last development [80]. Thus, at present it is 
not possible to distinguish between whether a direct interaction of Mg-Proto IX or heme activates 
transcriptional NGE, and whether Mg-Proto IX/heme only form intermediate molecules in the 
retrograde signal transductbn pathways. In contrast, 3'-phosphoadenosine 5'-phosphate (PAP) levels 
accumulate in response to drought and high light stresses and have been shown to alter RNA 
metabolism by directly inhibiting 5' to 3' exoribonuclease activity, leading to increased NGE [73]. 
The isoprenoid precursor MEcPP has also been identified as a retrograde signaling molecule 
through identificatbn of a genetic screen designed to identify genes involved in the regulation of 
HYDROXYPEROXIDE LYASE (HPL), a nuclear stress inducible gene that encodes a plastidial 
bcalized protein leading to adaptatbn responses upon bbtic stress conditions [74]. The resulting cehl 
mutant encodes for l-HYDROXY-2-METHYL-2-(E)-BUTENYL-4-DIPHOSPHATE SYNTHASE 
(HDS), responsible for catalyzing the conversion of MEcPP to hydro xymet hy lb utenyl diphosphate 
(HMBPP) in the plastidial methylerythitol phosphate (MEP) pathway. In contrast to the gun mutants, 
cehl does not coordinate PhANGs expression. Rather, the accumulation of MecPP increase salicylic 
acid (SA) levels through modulation of the SA biosynthetic ICSl transcript abundance and lead to 
Pseudomonas resistance [74]. 
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Thus, while a number of putative metaboUc mediators of plastidial retrograde signaling have been 
identified to date (based either directly on metabolite accumulation patterns or through genetic screens 
identifying biosynthetic coirponents), technical advances and approaches to study and elucidate the 
precise molecular mechanisms are still largely limited [71]. Furthermore, multiple independent 
signaling pathways might operate in parallel, or are only switched on during specific perturbations 
which increases the complexity associated with understanding and integrating or^nellar crosstalk. 
Thus, while the identification of metaboUc changes is possibly the first steps towards elucidating 
certain chemical or molecular signaling pathways that might be involved in retrograde signaling, more 
rigorous analytical approaches and subcellular advances [71 ] will provide the validation and context of 
metaboHc changes associated with organeUar dysflinctba 

3.3. Sugar Signaling 

Not all metaboHtes have a direct infiuence on the molecular network; alterations in the endogenous 
metaboUte levels could also reflect a change in the molecular flmctbn of its gene products. Clearly, the 
best and most extensively studied example of this is the sugar signaling network. Sugars particpate in 
glycolysis as the main respiratory carbon source to fliel the energy svpply of the cell. The first 
committed step towards mitochondrial respiratbn involves the phosphorylatbn of glucose and fructose 
to their respective hexose phosphorylated C6 intermediates via HEXOKINASE 1 (HXKl) or isomerization 
of the aldose or ketose phosphorylated CI derivatives to their respective hexose-6-phosphate moieties. 
HXKl is a conserved orthologous protein of the outer mitochondrial membrane [81,82] linking cytosoHc 
glycolysis to mitochondrial metabolism [83,84]. A loss of ftinction screen identifying glucose-insensitive 
(gin) mutants has elucidated HXKl as an important regulatory node in sensing internal glucose levels [85]. 
In this regard, HXKl localizes to the nucleus and activates nuclear gene expression through protein 
complex formation with the vacuolar H^ ATPase Bl (VHA-Bl) and 19S regulatory particle of the 
proteasome subunit (RPT5B) [85]. The biosynthetic activity of HXKl can be uncoupled from its 
sensing role, as catalytic inactive variants can complement the glucose sensitive pheno types, including 
reduced shoot and root growth, delayed flowering and senescence, and altered sensitivities to the growth 
hormones auxin and cytokinin [86]. However, low glucose levels can also mediate sugar signaling 
pathways independent of HXKl signaling [87] (see also bebw). During a range of abbtic stresses (see 
Table 1 for specific stresses), endogenous sugar levels change significantly across a range of stress 
stimuli (Table 1 ). 
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Table 1. Environmental perturbations lead to alterations in soluble sugar levels, gene 
e^qjressbn patterns and enzyme activities mArabidopsis thaliana. 



Abiotic Soluble 
stress sugar 



Sugar 
response 



Gene 
expression 



Gene 



Gene 
model 



Enzyme 
activity 



Enzyme Reference 



Cold 



glue, 
fruc 



increase 



increase 



FXK3 
GOLS3 
HXK2 
RS5 

SPSl 
SPSIF 



At5g51830 
Atlg09350 
At2gl9860 
At5g40390 
AtSgllllO 
At5g20280 



nd 



nd 



decrease SSR Atlg79440 



[88] 



increase 



Nitrogen 

(N) 
deficiency 



Xsugar 
(sue + 
glue + 
fioic) 



TPPB 
PGM 
A/N-Inv 
D 

G6PDH 
PGDH 
ABM 



Atlg78090 
Atlg78050 

Atlg22650 

Atlg24280 
Atlg64190 
At2g40220 



increase 



GlueoK 
FUM 



increase 



decrease 



TPS3 
CWINV4 



AtlglVOOO 
At2g36190 



decrease 



AGPase 
NADP- 
GAPDH 
NADP- 

IDH 

NR 

GS 
AspAT 



[89] 
[90] 



Potassium sue, 
(K) glue, 
deficiency fi^lc 



increase 



nd 



nd 



increase 



FrucK 
NADP- 

ME 
GS 
GOGAT 



decrease 



NR 



[91] 







PPCKl 


Atlg08650 






PPCK2 


At3g04530 






BAMS 


At4gl5210 






GWD3 


At4g24450 






GBSl 


Atlg32900 


Phosphate 


sue, 


increase GPT2 


Atlg61800 


(P) 


glue. 


increase ADG2 


At5gl9220 


deficiency 


fioic 


APL3 


At4g39210 






SPS4 


At4gl0120 






SPPl 


Atlg51420 






SUS3 


At3g43190 



increase 



cFBPase 

SPS 
PEP Case 
GS 
GDH 
ShiDH 



decrease 



FLNl 
FLN2 



At3g54090 
Atlg69200 



decrease PFP 



[90] 
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Abiotic 
stress 


Soluble Sugar 
si^ar response 


Gene 
expression 


Gene 


Gene 
model 


Enzyme 
activity 


Enzyme 


Reference 


Carbon 
availability 


Increase 
glue (exogenous 
feeding) 


increase 


ABM 
ABI5 
CTRl 


At2g40220 
At2g36270 
At5g03730 


nd 


nd 


[92] 







AMYl 


At4g25000 












BAMS 


At4gl5210 










sue, 
tre, 

increase 

glue, 
fruc 


C/VIFl 


Atlg47960 








Combined 
heat and 


FLNl 
decrease G6PDH6 


At3g54090 
At5g40760 


nd 


nd 


[93] 


drought 


HXK2 


At2gl9860 










PSL5 
SPSIF 
SPSl 


At5g63840 

At5g20830 
At5g20280 

















NDBl 


At4g28220 










sue. 






VHA-A 


Atlg78900 








Heat 


glue, 
fruc 


increase 


increase 


GOLSl 
A/N- 
INVC 


At2g47180 
At3g06500 


nd 


nd 


[93] 




sue, 
glue, 
fruc 






APL3 


At4g39210 








Drought 


increase 


increase 


SPSl 
SUSl 


AtSgllllO 
At5g20830 


nd 


nd 


[93] 








G6PDH6 


At5g40760 
















CWINVl 


At3gl3790 
















CINVl 


Atlg35580 














increase 


GDHl 


At5gl8170 








Drought 


sue 


increase 


GDH2 

VAC- 
INV 


At5g07440 
Atlgl2240 


nd 


nd 


[94,95] 








decrease 


ADKl 
GDH3 


At3g09820 
At3g03910 








Anoxia 


glue, 
fruc 


increase 


increase 


ADHl 
SUS4 


Atlg77120 
At3g43190 


nd 


nd 


[16] 
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Table 1. Cont. 



Abiotic Soluble Sugar Gene 
stress sugar response expression 



Gene 



Gene 
model 



Enzyme 
activity 



Enzyme Reference 



increase 



Salinity 



sue, 
glue 



increase 



PYD4 
BAMl 
GPAT5 



At3g08860 
At4gl5210 
At3gll430 



nd 



nd 



[96] 



decrease 



PDF 1.2b At2g26020 
PDF 1.2 At5g44420 



Abbreviatbns: ABI; ABA insensitive, ADG; ADP glucose pyrqjhospliorylase, ADH; alcoliol deliydrogenase, AMY; alpha- 
amylase, A/N-ESfVC; alkalineAicutral invatase C, APL; ADP -glucose pyrophosphoiylase, BAM; beta-anylase, cFBPase; 
cytosolic fructose bisphosphatase, CATF; cell wall/vacuolar inhibitor of fructosidase, CWTNV; cell wall invertase, CDSfV; 
cytosolic invertase, SSR; succinic semialdehyde reductase, fruc; fructose, FrucK; fructokinase, FXK; fructokinase, FLN; 
fructokinase-like, FUM; fumarase, ^uc; glucose, GBS; granule-bound starch synthase, GlucoK; gluconokinase, GOGAT; 
ferredoxin-glutamine-2-oxoglutaratc aminotransferase, GOLS; galactinol synthase, GPT; glucose 6 phosphate/t)hosphate 
transporter, OS; glutamine synthetase, GWD; glucan-water dikinase; G6PDH; glucose-6-phosphate dehydrogenase, 
GDH; glutamate dehydrogenase, GOLS; galactinol synthase, GPAT; glycerol-3-phosphate sn-2-acyltransferase, HXK; 
hexokinase, NADP-ME; NADP-dqjendent malic enzyme, nd; not determined, NDBl; NADH dehydrogenase Bl, NR; nitrate 
reductase, PDF1.2; plant defensin 12, PDF1.2b; plant defensin 1.2B, PEPCase; phosphoenolpyruvate carboj^flase, PFP; 
pytcp hosp hate- fructose-6-phosp hate phosphotransferase, PGDH; 6-phosphogluconate dehydrogpnasej PGM; phospho^ucortutase, 
PPQC; p hosphoeno|5yruvate carbojylase kinase, PSL5; priority in sweet life 5, PYD4; pyrimidine 4, RS5; rafiinose synthase 5, 
ShiDH; shikimate dehydro^nase, SPP; sucrose phosphate phosphatase, SPS; sucrose phosphate synthase, SPSIF; sucrose 
phosphate synthase IF, sue; suaose^ SUS; sucrose synthase, tre; trehalose, TPP; trehalose-6-phosphate phosphatase, TPS; 
thrdialose-6-phosphate synthase, VAC-INV; vacuolar acid invertase, VHA-A; vacuolar ATP synthase subunit A. 



Prevailing sucrose levels are sensed by the SUCROSE NON FERMENTATION 1 RELATED 
PROTEIN KINASE 1 (SnRKl), a conserved heterotrimeric protein kinase protein complex consisting 
of KINIO and KINll subunits. KINIO is the orthologue of the mammahan AMP activated protein 
kinase (AMPK), while KINll is the orthologue of the yeast sucrose non fermenting (SNFl) protein 
kinase, activated in response to low cellular glucose levels [97]. Plant SnRKls modulate sugar and 
ABA signaling pathways in response to energy depleting perturbations by controlling stress responsive 
gene expression [98]. Thios, while KINIO and KINll are predominantly localized in the cytoplasm, upon 
carbon and energy limitation, and form a protein complex that is sequestered to the nucleus where it has 
key transcriptional regulating capabilities. SnRKl transcriptional activity is inhibited by the levels of the 
endogenous phosphorylated sugars, trehalose-6-phosphate (T6P) and glucose 6-phosphate (G6P) [99-101], 
and upon phosphate starwtbn [102]. The phosphorylated sugar, T6P, also promotes thioredoxin mediated 
redox activatbn of ADP-glucose pyrophosphorylase (AGPase) involved in starch bbsynthesis in response 
to cytosolic sugar levels [103]. Overexpressbn of the basb leucine z^yper domain transcription fector, 
bZIPll further leads to the accumulatbn of endogenous T6P levels [104]. However, sucrose also 
represses the transktbn of bZIPl 1 leading to sucrose regulated changes in amino acid metabolism [105]. 
Multple lines of evidence now exist for the direct correlation between T6P and sucrose levels (e.g. [105,106]), 
and this association has been greatly accelerated by the accurate identification and quantification of the 
metabolites in questioa T6P also acts as a shoot apbal meristem (SAM) localized sugar signal for 
regulating the onset of flowering in response to environmental perturbations, including day length, 
teirperature, hormonal cues, and carbohydrate availability [106]. KINIO overexpression fijrther exhibit 
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glucose hypersensitivity (see also below) and alters carbon and nitrogen metabolismby controlling the 
expressbn of certain genes such as nitrate reductase (NR) and AGPase [98]. Upon hypoxic conditions, 
heterobgous KINIO overexpression leads to increased bngevity of rice and Arabidopsis under 
submergence due to increased expression of afcohol dehydrogenase (ADHl) and pyruvate decarboxylase 
(PDCl) [107]. The expression of these genes can be restored by sucrose addition [107]. 

Interestingly, KINIO also interacts with CYCLIN DEPENDENT KINASE E; 1 (CDKE;1), a 
regulator identified in the mitochondrial ALTERNATIVE OXIDASE (AOX) dependent retrograde 
signaling pathway and part of a nuclear transcriptional mediator coirplex [108]. Orthobgues of both 
KINIO and other sub units of the mediator conplex have also been identified inNeurospora crassa to 
regulate AOX expression [109], suggesting an evolutionary conserved response in crosstalk between 
components of retrograde and sugar signaling (see also below). While the exact molecular mechanism 
is not clear, cdke;l mutants display increased sugar levels upon unstressed growth conditions, while 
aoxla and kinlO mutants have reduced sugar levels compared to their respective controls diiring the 
photo synthetically active period (Radomiljac and Whelan, unpublished results). The aoxla, kinlO and 
cdke;l mutants all show normal developmental phenotypes grown under non stressed conditions, 
suggesting that the signal cascade opemting between sugar signaling, retrograde signaling, and phenotypical 
alterations upon abiotic stresses can be elucidated in reverse genetic stiodies combined withphenomic 
and metabolomics approaches. While this has not been addressed to date, it will be an exciting prospect 
in future in order to use metabolomics directed ^proaches to understand mitochondrial retrograde signaling 
and metaboHc activities to enhance plant fitness and survival during limiting growth conditions. Furthermore, 
ABSCISIC ACID INSENSITIVE 4 (ABM) has also been identified as a central regulator of both 
mitochondrial AOX dependent [1 10] and chloroplast dependent [1 1 1] retrograde signaling, modulating 
AOXla responses during high Ught and drought stress [112]. An extensive crosstalk between ABM 
and sugar signaling has been demonstrated [92], while abi4 mutants display increased salt tolerance [113]. 
Thus, while organellar retrograde signaling share a significant degree of crosstalk [114], metabolic 
changes related to sugar metaboHsm, signaling, and responsive gene expression further impHcates a 
higher degree and extent of intracellular small molecule communication during abiotic stresses. In 
additbnto the influence of sucrose and glucose mediated signaling, fructose levels can also be sensed 
in planta. Screening several enzymes in fructose metabolismby luciferase reporter screens in transient 
Arabidopsis mesophyll protoplasts have identified fructose insensitive (fin) mutants based on their 
inability to affect photosynthetic gene expression [115]. finl encodes for the metabolic enzyme 
FRUCTOSE- 1,6-BIPHOSPHATASE (FBPase) [115], responsfcle for fructose- l,6-b5)hosp hate steady 
state levels. Inhibition of cytosolic and plastidial FBPase activities severely limit photosynthesis in 
potato plants [116,117]. In addition, plastidial focalized FRUCTOKINASE-LIKE (FLN) proteins are 
targets for a plastidial thioredoxin, TRX z, that regulate plastid-encoded RNA polymerase (PEP) 
dependent transcription in plastids [118], suggesting that subcellular sugar metabolism and/or signaling 
have profound effects on coordinated nuclear organellar responses. 

TARGET OF RAPAMYCIN (TOR) proteins are evolutionarily conserved Ser/Thr kinases that 
control growth signaling pathways in response to stress signals through moduktbn of ribosome biogenesis, 
translation and primary metabolism [119]. Upon osmotic stress, TOR interacts with RAPTORl 
(a regulatory protein of TOR) to regulate the activity of S6 kinase, a protein conferring osmotic 
hypersensitivity to transgenic Arabidopsis plant [120]. TOR also interacts with LST8 proteins that 
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modulate plant growth, flowering, nitrate assimilation and sugar metabolism [121]. Although, the most 
recent player to be involved in integrating cellular sugar and energy depleting statuses, TOR signaling, 
in combination with SnRKl signaling, has been implicated as an important regulatory component 
during mitochondrial retrograde signaling. As kinases, both SnRKl and TOR proteins are involved in 
the post-translational phosphoryktbn modification of other target proteins. The alteratbn in the 
phosphorylatbn status due to perturbations in mitochondrial membrane potential has also been proposed 
to be a mechanism for the transbcatbn of putative transcription fectors encountered on the outer 
mitochondrial membrane [81,82] to the nucleus [122]. In addition to the identification of the kinase 
CDKE;1 as a regulator of AOXl a e5q)ression (but exhibiting post- transcriptional regulation of AOXla 
protein abundance) [108], an enigmatic questbn is the identification of the main metabolic signals that 
activate kinases and transduce the phosphorylatbn status. Furthermore, the overlap between sugar signal 
kinases and the retrograde kinases are a striking convergence point that needs to be explored in future. 

Thus, while sugar accumulatbn could potentially increase the osmotic potential within the cell 
(see introduction) or fLiel respiration through glycolysis, the TCA cycle and mitochondrial oxidative 
phosphorylation, a third flindamental consequence is the participation of sugars in specific signal 
transduction mechanisms during abiotic stress conditions. The specificity of the sug&r signaling 
response (ie. glucose, Iructose and sucrose signaling can be distinguished from each other) attests to a 
functional role for the individual sugars apart from an osmoprotectant or respfratory intermediate. In 
regards to the crosstalk and mechanisms operating between the different sugars, and the interactions 
within the rest of the cell; these are exciting research avenues where metabolomics driven endeavors 
can greatly provide novel insight into the molecular mechanisms underlying the specific abiotic stress 
responses. A key questbn arising fi-om the data presented in Table I (and in context with this review) 
is the estabHshment of whether the same core su^r signaling pathways are activated and/or transduced 
during the different abiotic stresses. Despite the individual sugars showing a striking similarity in 
response to a wide range of stresses, with multiple marker genes correlating with these (Table I), 
addressing the robustness in the response of sugar signal perturbations (and the vpstream cascade of 
events) is still largely lacking. The parallel time course analyses of the perceptbn and initiatbn of 
su^r signaling, as well as the molecular players activated vpon specific abiotic stresses, will he^ 
elucidate some of the key coirponents of this signaling response. 

One other aspect to address during the evaluation of metabolite changes is the robustness or 
diversity in metabolic activities and steady state metabolite abundances on species and developmental 
level. While the aforementioned disciission has exclusively focused on the molecular pathways or 
mechanisms associated with metabolite linked changes in the model species, Arabidopsis thaliana, a 
great deal of literature has also accumulated in other important crop and forage species where 
expansive metabolomic experiments have been performed upon various specific abiotic stresses (for 
example, salinity and drought stress [123-127]). This data could help to evaluate the trans latability of 
metabolic changes observed, and more importantly, the associated molecular mechanism(s) underlying 
the response, from model to inportant crop species. As an example, combined ionomics, transcriptomic 
and metabolomics studies on Lotus genotypes illustrated that the identification of genotype-specific 
molecular responses may provide a more reaHstic ratbnal for the design and implementation of 
tolerance traits under salt stress. The differential display of osmol)^ and other small molecule 
accumulation vpon drought stress of different Lotus genotypes suggest that coordinating metabolic 
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activities are orchestrated in a finely tuned manner [125,126], and might need more global approaches 
to understand these adaptive mechanisms (see also below). Similarly, investigating metabolite changes 
in different maize hybrids in greenhouse grown conditbns do not correlate with water deficit tolerance 
traits observed under field conditions [128], suggesting that finding distinct metabolic cues and under tying 
activities require molecular biology and agronomical approaches to converge in order to find practical 
field based solutions to improve plant growth and performance vpon non optimal growth conditions. 

4. Integration with the Other — Omics Platforms: Phenomics, Genomics, Proteomics, 
and Metabolomics 

While this review has largely focused on transcrpt- metabolite cohort changes, and the signal transduction 
pathways they participate in, metabolomics can be applied to a diverse range of -omic platforms to 
gain flinctbnal insight into physiobgical and molecular responses during abiotic stresses. The 
application of metabobmics to phenomics, genome wide associatbn (GWA) and recombinant mapped 
populatbns, genome and quantitative trait loci (QTL) sequencing efforts have provided the timely 
opportunity to explore the molecular and bbchemical basis and pathway interactions for coirp lex plant 
abbtic stress interactbns [48,129]. As an example, in a GWA mapping strategy appHed to 289 maize 
inbred Unes, the bbchembal conpositbn of 26 metabolites were strongly associated with single nucleotide 
poljmiorphisms (SNPs), including 15 distinct SNP-metabolite correlations that underpinned the genetic 
variance, observed in these hybrids [130]. Such discovery modules will greatly accelerate the fiiture 
directives to understand the complexity associated with, and strategies to improve, plant performance 
upon specific abiotic stresses. Fov Arab idopsis accessions, a similar approach has been applied upon C 
and N limitation and, once again, suggests that the resulting metabolic networks are highly variable 
and specific in terms of the limiting resource [131]. Similarly, upon cold stress, ecotype variance in 
Arabidopsis accessions could identify specific lipids or lipid classes (e.g. glucose/galactosylceramide 
(dl8:l/c24:0), monogalactosyldiacylglycerol (MGDG) 34:2, 34:3, 36:7), and triacylglyceride (TAG 
52:0 and 58:9)) through ultrahigh performance liquid chromatography Fourier transform MS (UPLC FT 
MS) technology positively correlated with low teirperature acclimatbn; albwing for potential metabolic 
bbmarker discovery that can be used in breeding strategies [132]. 

Furthermore, while perhaps less suited to proteomic studies dealing primarily with relative protein 
abundances, the integration of metabobmics with quantitative proteomics and, especially, enzyme 
activities [48], have also yielded fer greater insight into the (post-) trans lational metaboHc regulation of 
plants vpon, for exairple, limited carbon availability [48]. MetaboHtes can serve as transient 
modulators of protein fijnctbn and therefore have useful applications in the dynamic signal 
transduction pathway that lead to the activation of abiotic stress tolerance [133]. Cross et al. suggests 
that altered biomass and growth limitations of 24 genetically diverse Arabidopsis accessions ipon 
carbon deprivatbn is not a result of absolute metabolite levels of starch, sugars, or amino acids; but, based 
on the measured enzymatic capacity, related to the relative flux through the metabolic system [48]. As 
such, fluxomics, currently defined as the measurement of isotopic (stable or radioactive) label enrichment 
of metabolites (or pools) across time in organs, tissue or cell types, plays a pivotal role in linking 
changes in metabolite levels to that of genomic and proteomic platforms. The technologies and 
approaches to measure these on a more comprehensive, high throughput basis are still being developed 
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(for recent conprehensive review, [134]). One major limitation is our largely fragmented knowledge 
about plant metabolism. The identification of "new" metabolites or metabolite derivatives (both in 
primary and specialized metaboHsm), as well as novel metabolic bbsjoithetic pathways, is inperative 
to interpret metabolic flux labeling and account for flux balance equatbns. Current flux methodologies 
rely thus heavily on eirpirical mathematical modeling to obtain a reasonable goodness of fit of 
obtained experimental isotopic distributbn to that of curated metabolic models. As a result, fluxomics 
are currently more speciflc and applicable for tailoring experiments to resolve specific fluxes with high 
precision by labeling with precursor isotopes that are in cbse proximity to the area or flux distribution 
point of interest, or within a well-defined pathway. MetaboHc flux can also be only truly reflective of 
in vivo conditions if each metaboHte (and its resulting label enrichment) in the pathway, starting with 
total label iptake and utilization, can be quantitated, and these parameters need to be readily measured 
and be accounted for in any mass calculation equation. One way to currently gain biological insight 
into dynamic metabolite states is to define the degree of label enrichment within each individual 
metabolite (as assessed via hyphenated MS technologies), without interfering flux dynamics from it. 

As an example, in Arabidopsis roots exposed to oxidative stress; '^C glucose isotopic labeling 
incorporation into TCA cycle organic acids, and sucrose was significantly decreased, although the absolute 
concentration of these respective metabolites remained unaltered [135]. Upon recovery, experiments 
from oxidative stress, the metabolic recovery shifts (both steady state and dynamic metabolite levels and 
enrichment, respective^) precedes the enzymatic activities of the majority of the biosynthetic gene 
products [136], suggesting that combined metabobmics/fluxomics and proteomic approaches act as 
important platforms that can pinpoint the metabolic step(s) involved inplant stress responses. 

Another major area of advancement of metabolomic data through the combination of fluxomics, is 
the broader appHcation of isotopes to plant systems grown in closed faciHties svppUed with different 
isotopic labeling precursors (^^C, ^"^S) [34,137-139]. This approach leads to a greater confidence 
in current metaboHte annotation and also promise to expand the plant metabolite depository needed to 
improve current flux balance models. Greatly aiding these assignments is inproved instrumentation, 
including Fourier transform ion cycbtron resonance mass spectrometry (FTICR MS) and Orbitrap 
mass spectrometers (see also Section 2), that albws for mass accuracy and isotope elemental coirpositbn 
determinations. Combining such larger scale metaboUte/flux studies have not been applied to plants ipon 
abiotic stress, however, strategies to incorporate this promise to uncover unique and diverse metaboUte 
identities and signatiires that participate in various signaling and metabolic roles. 

Linking chemical genomics, fluxomics and metabolite profiling technologies have the further 
potential to uncover novel signal transduction pathways, and also for genetic engineering strategies to 
transiently perturb transcript or protein levels in response to stress stimuli (reviewed in [24,140]). In 
addition to these experimental integration strategies, subcellular metabolomics [134,141], tissue specific 
metabolite profiling [142], and the spatial compartmentation of cell specific transcriptional changes upon, 
for example nitrogen and iron deficiency [143, 144], have also provided more technical advances in refining 
models and hypothesis surrounding microscale metaboUc changes upon nutritional stresses to date. 
Combined with forther technologies, including fluxombs and bnombs, multidisciplinary approaches to 
understand abbtb stress adaptation and tolerance are rapidly expanding, and the next challenge is the 
systematic integration of metaboUc data into its correct molecular context. 
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5. Conclusions 

Metabobmics has taught us a great deal about the diversity in and dynamics of systems driven 
bb logical approaches during abiotic stresses. Ranging from osmoprotectants to allosteric regulators of 
protein properties, as well as direct interactors and modulators of DNA and RNA structures to alter 
nuclear gene expression, small molecules and the study thereof (metabolomics) have illustrated that, in 
combination with other -omic datasets, metabolites, proteins, and transcripts do not always "play by the 
rules of the game". While it challenges m to shift our assumptions and interpretation of metabolic changes 
in terms of how we think metabolism operates, it also greatly attests to the complexity within biological 
systems. Our current understanding of data structure relationships (particularly transcript-metabolite 
or metabolite- metabolite comparisons), and resulting hierarchical regulation structures, have been 
significantly improved by addressing metabolomics driven questions by careful designing experiments 
that will aid in understanding the complexity and regulation mechanisms during abiotic stresses. Plant 
metabolomics has taught us that, while metabolomics is a powerful tool to understand specific abiotic 
stresses (and even how these stresses overlap), we need the integration of other platforms or genetic 
variability in order to understand and integrate from what we measure to what the relevant molecular 
context is of these metabolic changes. However, metabolomics as a standabne entity will still provide 
valuable information regarding the identification of uncharacterized metabolites and/or derivatives, 
greater quantification of metabolites on both tissue specific and subcellular levels, as well as the dynamic 
movement of soluble metabolites between cells and tissues, especial^ vpon a range of stresses. In conclusion, 
metabobmics is a crucial link to decipher the molecular fijnctional mechanisms(s) underlying specific 
abbtic stress responses. 
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